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Abstract 

Based on the (3+l)-dimensional hydrodynamic model, the space- 
time evolution of hot and dense nuclear matter produced in non- 
central relativistic heavy- ion collisions is discussed. The elliptic flow 
parameter V2 is obtained by Fourier analysis of the azimuthal distri- 
bution of pions and protons which are emitted from the freeze-out 
hypersurface. As a function of rapidity, the pion and proton elliptic 
flow parameters both have a peak at midrapidity. 
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One of the main goals in relativistic heavy-ion physics is the cre- 
ation of a quark-gluon plasma (QGP) and the determination of its 
equation of state (EoS) It is therefore very important to study col- 
lective flow in non-central collisions, such as directed or elliptic flow 
Q. Recently experimental data concerning collective flow in semi- 
central collisions at SPS energies has been reported ^. This 
data should be analysed using various models. Some groups 0, ^ 
have used their microscopic transport models to analyse the collective 
flow obtained by the NA49 Collaboration j^. In this paper we inves- 
tigate collective flow, especially elliptic flow, in terms of a relativistic 
hydrodynamic model. 

In non-central collisions elliptic flow arises due to the fact that the 
spatial overlap region of two colliding nuclei in the transverse plane has 
an "almond shape". That is, the hydrodynamical flow becomes larger 
along the short axis than along the long axis because the pressure 
gradient is larger in that direction. Therefore this spatial anisotropy 
causes the nuclear matter to also have momentum anisotropy. Con- 
sequently, the azimuthal distribution may carry information about 
the pressure of the nuclear matter produced in the early stage of the 
heavy-ion collisions p. 

The relativistic hydrodynamical equations for a perfect fluid rep- 
resent energy-momentum conservation 

S^T'^" = 0, (1) 
T'"' = {E + P)u^v!' -Pg^"" (2) 

and baryon density conservation 

d^n^ = 0, (3) 
< = nen^, (4) 

where E, P, n-Q and are, respectively, the energy density, pres- 
sure, baryon density and local four velocity. We numerically solve 
these equations without assuming cylindrical symmetry [|^, ^] by 
specifying the model EoS and we obtain the space-time dependent 
thermodynamical variables and the four velocity. 

We use the following models of the EoS with a phase transition. 



Hagedorn's statistical bootstrap model |12| with Hagedorn tempera- 
ture Th = 155 MeV is employed for the hadronic phase. We directly 
use the integral representation of the solution of the bootstrap equa- 
tion [^] instead of using the very famous hadronic mass spectrum, 
exp(m/TH), which is the asymptotic solution of this equation. It is 
well known that this model has a limited temperature range, i.e., the 
energy density and pressure diverge at Th. This singularity, however, 
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disappears when an exclude volume approximation |]T4[ (with a Bag 
constant = 230 MeV) is associated with the Hagedorn model. In 
the QGP phase, we use massless free u, d and s-quarks and the gluon 
gas model for simplicity. The two equations of state are matched by 
imposing Gibbs' condition for phase equilibrium. Consequently we 
obtain a first order phase transition model which has a critical tem- 
perature Tq = 159 MeV and a mixed phase pressure of Pmix = 70.9 
MeV/fm^ at zero baryon density. 

We mention our numerical algorithm for the relativistic hydrody- 
namic model. It is known that the Piecewise Parabolic Method (PPM) 
1 15] is very robust scheme for the non-relativistic gas equation with a 



shock wave. We have extended the PPM scheme of Eulerian hydro- 
dynamics to the relativistic hydrodynamical equation. Note that this 



is a higher order extension of the piecewise linear method |16]. 

Assuming non-central Pb+Pb collisions at SPS energy, we choose 
very simple formulas for the initial condition at the initial (or passage) 
time to = 2ro/(7w) ~ 1.4 fm (rg, 7 and v are, respectively, the nuclear 
radius, Lorentz factor and the velocity of a spectator in the center of 
mass system) 

E{x,y,z) = Ei{z)e{zo-z)e{z + zo)p{rp)p{rt), (5) 
nB{x,y,z) = nBi{z)9{zo - z)e{z + zo)p{rp)p{rt), (6) 
v^{x,y,z) = i;otanh(2:/2;o) 

X 9{~zo- z)e{z+~zo)p{r^)p{rt), (7) 

where 6{z) is the step function, p{r) is the Woods-Saxon parameteri- 
zation in the transverse direction, 

P{r) = , (8) 

exp (r-m) + 1 

Ei[z) is Bjorken's solution jlTjj and the z dependence of the baryon 



density n-Bi{z) is taken from Ref. |18] 



Ei{z) = i?o X -^^4 , (9) 




't — z^ 

nBi(z) = K X 0.17-^^^^ . (10) 

See also Fig. |^. We have employed Bjorken's longitudinal solution 
just as an initial condition. This is in contrast to Ref. [^, |l^, in 
which Bjorken's boost-invariant solution was used as an assumption 
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Reaction Plane 



Transverse Plane 



Figure 1: Schematic view of the initial geometry in the center of mass system. 
The left figure shows the reaction plane and the right the transverse plane. 
The initial condition is in the region with slanting lines, b is the impact 
parameter vector, rp and rt are respectively the distances from the center of 
the projectile and the target nucleus in the transverse plane. 



and the hydrodynamical equation was numerically solved only in the 
transverse plane. 

At relativistic energies the Lorcntz-contracted spectators leave the 
interaction region after ~ 1 fm, we therefore assume the hydrody- 
namical description is valid only in the overlap region and neglect the 
interaction between the spectators and the fluid. Therefore we can 
say that our model gives a good description only in the vicinity of the 
midrapidity region and fails to reproduce directed flow at present. It 
may be possible to treat this problem if we use a hadronic cascade 
model for both spectators and particles emitted from the freeze-out 
hypersurface, together with the hydrodynamic model. 

There are four initial (and adjustable) parameters in our hydrody- 
namic model: 1) the energy density at z = 0, Eq = 2500 MeV/fm^, 
2) the factor in the baryon density distribution k = 2.5, 3) the ini- 
tial longitudinal factor e = 0.9 and 4) the "diffuseness parameter" 
5j: = 0.3 fm. In the present analysis we select these values 'by hand', 
i.e., we guess them. These parameters, however, should be chosen so 
as to reproduce the experimental data for the (pseudo-) rapidity and 
the transverse momentum distribution. To make our analysis more 
quantitative, we need this experimental data. We would like the ex- 
perimental group to analyze the centrality dependence of the hadron 
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Figure 2: Time evolution of pressure and baryon flow in the transverse 
plane. Left: The pressure contours. Right: The baryon flow velocity vector 
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spectra, especially, the (pseudo-)rapidity distribution. For this reason 
we wish to emphasize that our numerical results presented below are 
only preliminary. 

Figure |2| shows our numerical results for the temporal behavior 
of the pressure (left column) and the baryonic flow (right column) at 
z = in the non-central Pb+Pb collision with impact parameter 6=7 
fm at SPS energy. Initially almost all matter in this plane is in the 
QGP phase and there is no transverse flow anywhere by definition. 
At t = to + 0.5 fm we see the shell structure corresponding to the 
mixed phase with the same pressure ~ 70 MeV/fm^, and the initial 
pressure gradient gives the baryons transverse flow. The QGP phase 
disappears at i = io + 1-0 fm and after that the mixed phase occupies 
the central region. There is still no transverse flow near the origin due 
to the absence of a pressure gradient. At about t = to + 5.0 fm all 
the nuclear matter initially in the QGP phase has gone through the 
phase transition and is in the hadronic phase. We can see from these 
figures that the shape of the nuclear matter is changing from almond 
(top figure on page 5) to round (bottom figure on page 6), and the 
elliptic flow reduces the initial geometric deformation. 

The numerical results of the hydrodynamical simulation give us the 
momentum distribution through the Cooper- Frye formula |20] with 



freeze-out temperature Tf = 140 MeV. The elliptic flow parameter V2, 
as a function of rapidity y, is obtained from the momentum distribu- 
tion 



My) 




dp 



(11) 



Before calculating V2 in non-central collisions with impact param- 
eter b = 7 fm, we checked the numerical error in our hydrodynamic 
model in central collisions. Since there is no special direction in the 
transverse plane for head-on collisions, ideally the elliptic flow vanishes 
in the infinite particle limit. Performing the numerical simulation with 
6 = fm, we obtain the value of V2 as less than 10~^ percent, therefore 
we can safely neglect the numerical error. Note that the numerical er- 
ror in the energy and baryon density conservation of the fluid is less 
than one percent in our analysis. 

Figure |^ shows our results for the rapidity dependence of elliptic 
flow for pions in different transverse momentum regions. These results 
show that elliptic flow rises with transverse momentum pt and 
has a peak at midrapidity. This seems to be in contrast with the 
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Figure 3: Rapidity dependence of elliptic flow for pion. Four curves cor- 
respond to the different transverse momentum regions. The midrapidity is 
2.92. 




Figure 4: Rapidity dependence of elliptic flow for proton. Three curves 
correspond to the different transverse momentum regions. Note that the 
integral region of transverse momentum is larger than for pions. 
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experimental data obtained by the NA49 Collaboration Q. Their data 
appears to be slightly peaked at medium-high rapidity. 

Our results for V2 for protons are shown in Fig. ^. We see the 
same behavior as for the pion case. We obtain a larger V2 for protons 
than for pions because we are integrating over a larger transverse 
momentum region. Since the initial parameters in our hydrodynamic 
model have been chosen by hand, we would like readers to not take 
these results quantitatively. 

In summary, we reported our preliminary analysis of elliptic flow 
in non-central heavy-ion collisions using the hydrodynamic model. 
We numerically simulated the hydrodynamic model without assuming 
cylindrical symmetry or Bjorken's boost-invariant solution, using the 
extended version of the Piecewise Parabolic Method which is known 
as a robust scheme for the non-relativistic gas equation with a shock 
wave. We presented the temporal behavior of high temperature and 
high density nuclear matter produced in Pb+Pb collisions with 6 = 7 
fm at SPS energy. Our preliminary results showed that the elliptic 
flow parameter V2 has a peak at midrapidity for both pions and pro- 
tons and increases with transverse momentum. Since there are some 
ambiguities in the initial parameters of our hydrodynamical model, we 
should fix these parameter using experimental data for the rapidity 
distribution in non-central collisions. If we regard the hydrodynami- 
cal model as a predictive one, we can choose initial parameters using 
results from a parton cascade model, such as VNI [22]. The study of 
these issues is a future work. 

The author is much indebted to Prof. I. Ohba, Prof. H. Nakazato, 
Dr. Y. Yamanaka and Prof. S. Muroya for their helpful comments, and 
to Dr. H. Nakamura, Dr. C. Nonaka and Dr. S. Nishimura for many 
interesting discussions. The numerical calculations were performed on 
workstations of the Waseda Univ. high-energy physics group. 
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